Levy flights as subordination process: first passage times 
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We obtain the first passage time density for a Levy flight random process from a subordination 
scheme. By this method, we infer the asymptotic behavior directly from the Brownian solution and 
the Sparre Andersen theorem, avoiding explicit reference to the fractional diffusion equation. Our 
results corroborate recent findings for Markovian Levy flights and generalize to broad waiting times. 



An important issue in the theory of stochastic pro- 
cesses is the problem of first passage pj H 0, : Its solu- 
tion is a key to the understanding of chemical reactions, 
stability of states of dynamical systems under external 
perturbations, extinction of populations, and many other 
problems in natural sciences. In normal Gaussian dif- 
fusion, a standard method of solution of the first pas- 
sage problems is the method of images due to Kelvin 
0- tJ However, for Levy flights (LFs) it has recently 
been demonstrated that the images method leads to a 
result that contradicts the Sparre Andersen theorem, ac- 
cording to which the first passage time density (FPTD) 
of a random walk process asymptotically follows the uni- 
versal f{t) ~ t~ 3 / 2 behavior for any symmetric distri- 
bution of jump lengths 0, 0, 0- This statement is of 
great importance since it not only points out possible in- 
applicability of Kelvin's method to superdiffusive jump 
processes, but also poses intricate questions of a correct 
continuous limit for Levy flights and their description 
within fractional diffusion or Fokker-Planck equations. 

LFs are a paradigm for anomalous stochastic processes 
with a wide range of applications such as chaotic dy- 
namics, processes in plasma, transport in micelles, or 
even quantum systems @ H E El El El El . They 
can be defined in terms of a random walk process with 
long-tailed jump length distribution X(x) ~ a tl /\x\ 1+ ' 1 
(0 < fi < 2) Alternatively, LFs can be considered 

as appearing from a Langevin equation with 5-correlated 

Levy noise El El El 

The characteristic function of an 

LF given by 



P(k, t) — exp 



(1) 



is of stretched Gaussian type, where P(k, t) = 
^{P(x,t)} = P(x,t) exp(ikx)dx is the Fourier 

transform of the probability density function (PDF) 
P(x, t) dBHEl. In what follows, the generalized diffu- 
sion constant K^' will be set to unity. In the limit fj, = 2, 
the characteristic function Q reduces to the Gaussian 
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P(k,t) = exp (— K^k 2 i) , the characteristic function of 
a standard random walk with Gaussian limit distribution 
and finite variance (x (t)> = 2K&H EH13, and finite 
higher order moments. For the general case < /i < 2, 
only fractional moments of the form 



exist with S < fi which can most easily be obtained by 
using properties of Fox ii-functions [2 l|. 

The fractional diffusion equation |lf lllll El E3 



dt d\x\" • 1 ' J 



governing an LF reflects the inherent long jumps _by_ the 
non-locality of the fractional Riesz-Weyl operator 

MM 



(3) 
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d^P{x,t) = 
d\x\» 



1 
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2 cos (tt/V2) T(2 - /i) dx 2 
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, \x-i\^ 
(4) 

for 1 < a < 2, and an analogous expression for < /i < 1; 
its Fourier transform &{d»P(x, t)/0|af|"} is -\k\^P{k,t) 
[Tel 122^. In the limit (j, = 2, the fractional operator Q 
reduces to the standard second-order partial derivative, 
as it should. 

In terms of the dynamical equation for the PDF, the 
formulation of the first passage problem in terms of an 
absorbing (Dirichlet) boundary condition requires a cut- 
off in the Riesz-Weyl integral (@J tJ. Contrasting the 
universal Sparre Anderson result for the FPTD, it was 
found that the probability density of first arrival at a site 
differs from the FPTD, and is explicitly dependent on the 
index /i of the LF 7]. In the present work, we present 
an alternative derivation of the FPTD for LFs based on 
a subordination to a regular random walk process of the 
kind developed in Refs. [HE!; leading us to a new, a 
priori unexpected twist: Although the PDF of the free 
process (without boundaries) is described by the frac- 
tional diffusion equation, Eq.JJJ), this equation does not 
uniquely describe the whole process, i.e. the exact prop- 
erties of its sample paths. Depending on how anomalous 
transport statistics are introduced, we will discuss the 
appropriateness or failure of the images method. 
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Let us consider a stochastic process subordinated to 
a discrete random walk with some operational time 
that is defined by a one-sided Levy law [2j|. The 
notion of subordination implies that the correspond- 
ing random process can be understood as follows [3|. 
The motion of the random walker can be parameter- 
ized by the number of steps n, the PDF of the random 
walker's position after n steps being given by the PDF 
Pnw(x,n). The corresponding characteristic function 
ip n (k) = J_ oc Prw(x, n) exp (ikx) dx is determined by 
the PDF of step lengths X(x) such that <p n {k) = [A(fc)]™ 
in terms of the characteristic function X(k) of the step 
length PDF, compare, for instance, Ref. Q- 

The number of steps n itself that may be considered as 
the internal operational time of the process, is a nondc- 
caying random function of the physical time (clock time) 
t. Denoting by p n (t) the probability to perform exactly 
n steps up to time t, we obtain 



P(x, t) = Prw(x, n)p n (t). 



(5) 



At long times we suppose that a continuous operational 
time T may be introduced instead of the discrete index 
n. Moreover, we assume that the continuous analog of 
Prw{x, n) exists. In this limit, we find from QSj 



P(x,t) 



P x (x,T)p T (T,t)dT, 



(6) 



where P x (x, T) is the PDF to be at x at operational time 
T and pt(T, t) is the PDF to be at operational time T at 
clock time t. This continuous limit corresponds exactly 
to the mathematical notion of the subordination, see Sec. 
X.7 of Ref. 0: 'If {X{T)} is a Markov process with con- 
tinuous transition probabilities and {T(t)} a process with 
non-negative independent increments, then {X(T(t))} is 
said to subordinate to {X(t)} using the operational time 

t: 

Instead of having the PDF Pnw(x,n) of a Brownian 
random walk, let us assume the more general case that 
Prw(x,ti) corresponds to an LF, so that Prw{%, "J, is 
given by a symmetric Levy distribution of the form 



P RW (x,n) = 



Ion 1 / 01 



,a,0 



(7) 



(0 < a < 2), where Iq is a scaling factor with length as 
physical unit. Note that the limit a = 2 corresponds to 
a Gaussian profile for Prw(x, n). Let us additionally as- 
sume that the number of steps per unit of the physical 
(clock) time t is distributed according to some distribu- 
tion with a power-law tail: p(n, At = 1) oc n~ 1 ~' 3 , with 
< (3 < 1. Then, according to the generalized central 
limit theorem, at long times p n (t) tends to a continuous 
limit distribution corresponding to the one-sided Levy 
law 



MP 



(8) 



Here, To is a scaling factor with physical unit time. The 
one-sided character of (JSJ) ensures that Pt(T < 0,t) = 

To obtain the limit distribution P(x, i) based on rela- 
tions 10) and J3J), we first Fourier-transform Eq. (0 in 
respect to x. With expression for the characteristic 
function of a symmetric Levy stable density, we find 



P(k,t) = J exp[~\klo\ a T)(Jf) 



i/0 



xL(V(^) 1//3 ;f3,-p)dT. (9) 

This exactly corresponds to the Laplace-transform of a 
one-sided Levy density with Laplace variable u = \k\ a lQ. 
The Laplace transform of a one-sided Levy stable density 
is known @,E1: PT( u ,t) — exp (— u^i/r ), and thus 



P(k,t) = exp (-\k\ al3 [Zq'V 1 ] t) . 



(10) 



This expression, in turn, is the characteristic function of a 
symmetric Levy stable density with index a(3. The scal- 
ing factor K^ a ^ = Iq 13 /tq can be interpreted as the asso- 
ciated fractional diffusion coefficient that depends on the 
indices of the corresponding subprocesses only through 
their product 7 = a(3. 

The PDF HI Up fulfills the fractional diffusion equation 
P|) with order /1 = 7. This can also be shown on gro unds 
of the subordination scheme developed in Refs. |23l l25j. 
However, it is remarkable that the dynamical equation 
© defines the PDF P(x,t), but not the process (or bet- 
ter, the identity of the subprocesses) itself; unless we have 
to do with the limiting case 7 = 2 that necessarily corre- 
sponds to a — 2 and [3 — 1, defining the process uniquely. 
In other words, although the PDFs of all processes with 
identical 7 are the same, these processes may still differ 
in the fractal dimension of their sample paths given by 
the set of jumps of random length corresponding to LFs 
with index a; and similarly they differ in the nature of 
the connection between operational time T and physical 
clock time t according to the one-sided Levy stable den- 
sity with index /?. This issue is of great importance for 
first passage problems, as we demonstrate now. 

Our subordination procedure corresponds to a (ran- 
dom) change of the time variable of the process from 
operational time T to clock time t. This allows us to 
solve a number of problems connected to the underlying 
random process without explicitly referring to the frac- 
tional equation ©. We here consider the first passage 
across a boundary located at x = 0. 

Let y n {xo) be the survival probability on the positive 
semi-axis (i.e. the probability of not crossing the bound- 
ary within the first n steps) after starting at xq > 
at n = 0. According to the Sparre Andersen theorem 
@ [J the asymptotic form of this probability does 
not depend on the jump length distribution if only it is 
symmetric. For large number of steps, one invariably has 
^ n (xo) — c(cco)"- _1 ^ 2 7 where the prefactor c(xq) depends 
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on the initial position xq, as well as on a and Iq. On sub- 
ordinating the number of steps n to physical clock time 
t we see that the survival probability up to time t is 



(11) 



where p n (t) is the probability that exactly n steps occur 
within clock time t. Note that equation (jl 1|) corresponds 
to the spatial average of Eq. J5J. Changing from n to 
the continuous operational time variable T, we get: 

/>oo 

y{t;x )= 9(T;x )p T (T,t)dT. (12) 
Jo 

This and \& n (xo) — c(a; )n -1 / 2 gives rise to 



y(t;xo) ~ c(xo) 



^(^) (?)""• M 



as derived in the Appendix. We now obtain the FPTD, 



/(<) 



dJ?(t;x ) 
dt 



c(*o) r d | M ^ 

tt 1 ^ ^ ^ 2/?y £1+1/(20) ■ 



(14) 

The following limiting cases can be distinguished: 

(i) If the subordination from operational time T to 
clock time t through px(T,t) is narrow with (3=1, i.e., 
p T (T,t) = 5(T-t/T ), the universal /(*) - t^ 2 behav- 
ior according to the Sparre Anderson theorem is recov- 
ered. In other words, in order to change this asymptotic 
behavior, one has to consider < (3 < 1 explicitly. 

(ii) If we consider the process subordinated to Gaussian 
diffusion [a = 2), then (3 = 7/2 and 



/(*) 



2c(z ) 

71-1/2^, 



r(i + 7) 



1/7 

£1+1/7 ' 



(15) 



This result has the same scaling as the FPTD derived 
through the method of images in Ref . • A few words on 
the interpretation of this seemingly paradox finding are in 
order. Result (|T5|l corresponds to a random walk process 
with a Gaussian jump length density X(x), so that the 
corresponding trajectory is that of a regular Brownian 
walk. It is therefore perfectly legitimate to use the images 
method for such a process, even though it is described 
by the fractional diffusion equation ©. In contrast, the 
genuine ('classical') LF discussed in Ref. Q corresponds 
to a broad X(x) ~ Zp/|a;|i +Q with a explicitly smaller 
than 2, but (3=1. For this strongly non-local process 
with fractal trajectory, the FPTD follows the result (i). 
In this latter case, that is, the method of images fails and 
the FPTD follows f(t) ~ t~ 3 / 2 . 

(iii) In general, for a given 7 one has to make sure that 
the inequalities 7/2 < (3 < 1 are fulfilled, since simul- 
taneously the two conditions (3 < 1 and (3 = j/a with 
< a < 2 have to be met. The FPTD for such a general 
7 therefore shows the following asymptotic behavior 



where < 7 < 2. In this scheme, the Sparre Andersen 
decay with exponent 3/2 is the slowest one possible. This 
makes perfect sense since due to the Levy stable form of 



Pn(t) 



-1-/3 



a broad distribution of single jump events 



f(t) <x r 



3/2 < S < 1 + I/7, 



(16) 



occurs in a finite time interval (0,t), increasing the like- 
lihood of crossing the boundary within any given finite 
time interval dramatically. The Levy nature of p n (t) thus 
leads to an oversampling of the space in comparison to a 
process (3 = 1. 

The subordination scheme developed here allows one to 
express the FPTD of a random process solely on the basis 
of the properties of the Sparre Andersen universality and 
the subordination map from operational time T to phys- 
ical clock time. Starting off from an LF in (n, t) coordi- 
nates, we introduce a broad distribution p n (t) ~ n~ 1 ~ 13 
of events per clock time interval At. We find that the 
resulting process is governed by the fractional diffusion 
equation Q3} whose order 7 = a(3 is a product of the Levy 
index of the jump length distribution, a, and the subor- 
dination distribution p n (t), (3. The knowledge of 7 alone 
is therefore insufficient to deduce the exact form of the 
jump length PDF and the trajectory it gives rise to. 

As a direct consequence, the resulting FPTD in the 
limit p n = 5 nt i (the trivial subordination with (3 = 1) 
fulfills the Sparre Anderson universality for any process 
with symmetric jump length distribution. This case in- 
cludes the case of genuine LFs as those discussed in 
Ref. 7], and is violated by the method of images for 
all < a < 2. Conversely, for a = 2, the subordination 
process has the trajectory of a normal Brownian random 
walk and is amenable to the images method to determine 
the FPTD. The result scales like /(*) - t- 1 - 1 /^ as previ- 
ously obtained. In the case of general a and (3, the range 
spanned by the exponent in the FPTD f(t) ~ t^ 1 ^ 8 is 
3/2 < 5 < 1 + I/7. 

We believe that above findings help in interpreting the 
inadequacy of the images method for genuine LFs as 
found in Q and moreover show that caution is neces- 
sary when generalizing well-known results from ordinary 
to anomalous diffusive processes. 

We acknowledge helpful discussions with Aleksei 
Chechkin and Yossi Klafter. 



Appendix A: CALCULATION OF THE 
( l/2)-ORDER MOMENT 

To obtain the final form of Eq. Q13[l . we have to eval- 
uate the integral 

1/9 /Tn\l/( 2 ' 3 ) / /T \ 1 /P 

which is equal to the (— l/2)-order moment 

m_ 1/2 {(3)= / r 1/2 m(3;-m 

Jo 

of the one-sided Levy stable density L(£;(3;— [3). To 
calculate this expression, one can follow two alternative 
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routes. As such techniques are useful when dealing with 
Levy processes, we sketch both ways. 

(i) The density L(£; (3; —0) is the inverse Fourier 
transform of the characteristic function y(fc) = 
exp (— Ifc^e - "' 3 / 2 ) for k > 0, continued to the negative 
semi-axis by tp(— k) = ip*(k). Then, 



M i -O'i = / ^Re 



00 r/fr 

— e - 4fc «exp 



-Z7T/3/2 



Interchanging both integrals, one can express the in- 
tegral over £ through Fresnel integrals to obtain 



-1/2, 

then — ifc to u we find for the overall integral: 



Changing 



M_ 



1/2 



09) 



"^Re 



dMexp(-u' 3 ))u- 1/2 . 



Noting that the integrand is analytical everywhere ex- 
cept for u = 0, vanishes for |it| — * oo in the lower right 
quadrant of the complex plane, and that the integral con- 
verges, we can change the contour of integration to the 
positive real axis, so that 



M_ 1/2 (/3) 



du exp I 



(Al) 



Substituting r\ = vP , we recover the the integral repre- 
sentation of the T-function, so that Eq. I|13l) follows. 

(ii) Alternatively, the density L(£; /3; —(3) is 
the inverse Laplace transform of exp (— u' 5 ) = 



n.i 



(0,1/0 
ion HI 



expressed in terms of the 



123 . By standard methods (26L 12 



Fox iJ-funct 

the inverse Laplace transform can be performed, yielding 



(0,1) 
(0,1/(3) 



After a substitution, the resulting integral 



M_ 1/2 (/3) 



-l/2rrl,0 
"1,1 



(0,1) 
(0,1/0) 



can be easily evaluated, by noting that it corresponds to 
the Mellin transform g(s) = §^°t s ~ 1 g(t) of H 1 \(z) at 
s = 1/2. The result can be directly identified with the 
definition of the ff-function [2lll2^ |. so that M_ 1/2 ((3) = 
/?- 1 r(l/[2/3])/r(l/2), reproducing Eq. JT^. 
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